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Abstract 

We analyze the recent experimental evidence for an excess of r-lepton production in several 
exclusive semileptonic S-meson decays in the context of two-Higgs-doublet models. These 
decay modes are sensitive to the exchange of charged scalars and constrain strongly their 
Yukawa interactions. While the usual Type-II scenario cannot accommodate the recent 
BaBar data, this is possible within more general models in which the charged-scalar 
couplings to up-type quarks are not as suppressed. Both the B — > D^*>tu t and the 
B — > tv t data can be fitted within the framework of the Aligned Two-Higgs-Doublet 
Model, but the resulting parameter ranges are in conflict with the constraints from leptonic 
charm decays. This could indicate a departure from the family universality of the Yukawa 
couplings, beyond their characteristic fermion mass dependence. We discuss several new 
observables that are sensitive to a hypothetical charged-scalar contribution, demonstrating 
that they are well suited to distinguish between different scenarios of new physics in the 
scalar sector, and also between this group and models with different Dirac structures; 
their experimental study would therefore shed light on the relevance of scalar exchanges 
in semileptonic b — > ct~9 t transitions. 



1 Introduction 



The BaBar collaboration has recently reported an excess of events in two semileptonic transi- 
tions of the type b — > ct~v t . More specifically, they have measured the ratios [l] 

Bt(B^Dt-V t ) BaBar avg 

R(D*) = Br ( B _^ D * T V A Ba J- 0.332 ±0.024 ±0.018 = g 0.354 ± 0.026 , (1) 
v ; Br(5 — >■ D*£~i?i) ' v ; 

which are normalized to the corresponding decays into light leptons £ = e, /i. The second value 
given in each line is the average with the previous measurements by the Belle collaboration (2j[3] , 
which also yield central values corresponding to an excess, but not significantly so. These 
relative rates can be predicted with a rather high accuracy, because many hadronic uncertainties 
cancel to a large extent. The Standard Model (SM) expectations [4jj7] are significantly lower 
than the BaBar measurements. If confirmed, this could signal new physics (NP) contributions 
violating lepton-flavour universality. 

A sizable deviation with respect to the SM prediction was previously observed in the leptonic 



decay B — > r u T , when combining the data by BaBar [8j|9] and Belle 10 . The world average 
Bt(B~ ->■ r _ P T ) = (1.65 ± 0.34) x 10" 4 |'xTJ used to be 2.5a higher than the SM prediction 
(0.733+o;o73) x 10~ 4 fl2) (taking the modulus of the Cabibbo-Kobayashi-Maskawa (CKM) 



matrix element \V u b\ from a global CKM fit). However, a value closer to the SM expectation 



13 



has 



been just reported by Belle 14 , leading to the new Belle combination Bi(B — > t z/^Beiie = 
(0.96 ± 0.22 ± 0.13) x 10 -4 , which we average with the combined BaBar result [9] to obtain 
Br(B~ -»■ t~v t ) = (1.15 ± 0.23) x 10- 4 . 

While more experimental studies are clearly needed, these measurements are intriguing 
enough to trigger the theoretical interest [4^|7| [l5}]20] . This kind of non- universal enhancement 
of the t production in semileptonic S-meson decays could be generated by NP contributions 
with couplings proportional to fermion masses. In particular, it could be associated with the 
exchange of a charged scalar within two-Higgs-doublet models (2HDMs), with the expected 
contribution to the transition amplitude being proportional to m T mb/M^ I± . This approach 
offers (obviously) a solution when considering scalar NP contributions model-independently; 
even in this general case, however, non-trivial predictions for other observables in these decays 
can be obtained. More specific models generally have difficulties in describing all the data. 
For example, the BaBar data on B — > Dt~v t and B — > D*r~i? T cannot be explained simul- 



taneously within the usually adopted Type-II scenario [1114,15 . We shall show, however, that 
they can be accommodated by the more general framework of the Aligned Two-Higgs-Doublet 
Model (A2HDM) [21], albeit creating a tension when including them in a global fit. 

Generic multi-Higgs-doublet models give rise to unwanted flavour-changing neutral cur- 
rent (FCNC) interactions through non-diagonal couplings of neutral scalars to fermions 
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The tree-level FCNCs can be eliminated by requiring the alignment in flavour space of the 



Yukawa matrices coupling to a given right-handed fermion 23 . This results in a very specific 
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structure, with all fermion-scalar interactions being proportional to the corresponding fermion 
masses, and implies an interesting hierarchy of FCNC effects, suppressing them in light-quark 
systems while allowing potentially relevant signals in heavy-quark transitions. The A2HDM 
leads to a rich and viable phenomenology 21 23 26 ; it constitutes a very general framework 
which includes, for particular values of its parameters, all previously considered 2HDMs with- 
out tree-level FCNCs |22|, and at the same time incorporates additional new sources of CP 
violation beyond the SM. 

In the following, we shall consider the phenomenology of b — > qr~i? T (q = u,c) transitions 
within a framework with additional scalar operators, assumed to be generated by the exchange 
of a charged scalar. Starting from the most general parametrization of such effects, we then 
specialize to various more specific models to examine their capability to describe the data and 
the possibility to distinguish between them. 

Our paper is organized as follows: In Sec. [2j we briefly describe the theoretical framework 
adopted in our analysis. In Sec. [3j we present our numerical results and show the parameter 
ranges needed to explain the present data. We proceed in Sec. [4]by analyzing various additional 
observables sensitive to scalar contributions, both integrated and differential, before concluding 
in Sec. [5] The appendices include a discussion of the relevant input parameters and details on 
the calculation for the semileptonic 5-meson decays. 



2 Theoretical framework 

We are going to assume that, in addition to the SM H^-exchange amplitude, the quark-level 
transitions b — > ql~v\ receive tree- level contributions from the exchange of a charged scalar. 
The effective low-energy Lagrangian describing these transitions takes then the form 



V2 

q=u,c l=e,fi,r 



C eS = v * E [h.Vm] + [q{gt l V L + gi l V R )b] [iVm]} , 

(2) 

where the effective couplings are, in the majority of 2HDMs, proportional to the fermion masses, 
g^idi ^ rriq u mi/Mfl±, g q ^ qdl ~ m qd mi/M^ I± . These explicit fermion mass factors imply negligi- 
ble effects in decays into light leptons (e, while decays involving the r receive potentially 
large contributions. Owing to the m q suppression, the coupling g q L does not play any relevant 
role in b — > ut~p t transitions, but it can give sizeable corrections to b — > ct~v t . 

We present next the key relations for including scalar NP contributions in the decays in 
question. When considering specific models, the main focus lies on the A2HDM, of which we 
give a short review. 



1 The obvious exception are leptonic meson decays, where the SM contribution is already suppressed by the 
light lepton masses, yielding a large relative contribution from the charged scalar. 
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2.1 b —± qr v T (q = u, c) decays 



Due to the helicity suppression of the SM amplitude, the leptonic decay B — >■ r u T is partic- 



ularly sensitive to the charged scalar exchange. The total decay width is given by 24,27 



r(B- 



— > T V T 



GWJl\V ub \ 2 



m B 



in 



mz 



m 



B 



:i + s c 



A T I 2 



(3) 



where S em denotes the electromagnetic radiative contributions, and the new-physics information 
is encoded in the correction^ 

( n qbl J M U2 
\9l - 9r ) m B <?= 



A 



9r 1 ml 



gb 



mi(m b + m q ) 



mi m b 



(4) 



absorbing in addition mass factors from the hadronic matrix elements. 

Semileptonic decays receive contributions from a charged scalar as well, but in this case the 
leading SM amplitude is not helicity suppressed; therefore, the relative influence is smaller. In 
addition, they involve momentum-dependent form factors. The B — > DWi decay amplitude is 
characterized by two form factors, f+(q 2 ) and fo(q 2 ), associated with the P-wave and S-wave 
projections of the crossed-channel matrix element (0|c7 M 6|5D). The scalar-exchange amplitude 



only contributes to the scalar form factor; it amounts to a multiplicative correction 24 



fo(q 2 ) fo(q 2 ) = fo(q 2 



tub — mr>) 2 



with 



(9l + 9 C R )( m B -m D ) 



cb 



(5) 



(6) 



mi[m b 

The sensitivity to the scalar contribution can only be achieved in semileptonic decays into 
heavier leptons. The decays involving light leptons can, therefore, be used to extract infor- 
mation on the vector form factor, reducing the necessary theory input to information on the 
scalar form factor. Furthermore, in NP searches, ratios of observables are usually used in order 
to make them independent of the normalization. This leads us to parametrize the different 
CP-conserving observables that we are going to consider in the form: 



O = c + CiRe^J + c 2 \5 T cbl 



(7) 



implying a discrete symmetry Im (5 T cb ) — > — Im(8l b ). The coefficients q, which contain the 
dependence on the strong-interaction dynamics, are parametrized in turn in terms of the vector 
form-factor slope p\ and the scalar density A(vr-vd)- For the former we use the value extracted 
from B — > Dtv {t = e, fi) decays. The function A(t> B • vp) oc fo(t)/ f+(t) has been studied on 
the lattice, in the range Vb • Vo — 1-1.2, and found to be consistent with a constant value 
A = 0.46 ±0.02, very close to its static-limit approximation (mg — mo) / (tub + tub) [7, 17,28 



This value is furthermore in agreement with QCD sum rule estimates 29 , 30 



2 Here and in the following we suppress in the notation the scale dependence of e.g. the quark masses and 
scalar couplings. 
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The decay B — > D*l~i?i has a much richer dynamical structure, due to the vector nature 
of the final D* meson. The differential decay distribution is described in terms of four helicity 
amplitudes, H±±, H 00 and H ot , where the first subindex denotes the D* helicity (±,0) and 



the second the lepton-pair helicity (±,0,t) in the .B-meson rest frame 31 32 . In addition 
to the three polarizations orthogonal to its total four- momentum q^, the leptonic system has 
a spin-zero time component (t) that is proportional to and can only contribute to the 
semileptonic decays for non-zero charged-lepton masses (it involves a positive helicity for the 
The corresponding Hot amplitude is the only one receiving contributions from the scalar 
exchange (H]: 

Hot(q 2 ) = l?g?V) (l-A^y (8) 



m 



The observables are then given in an expansion analogous to Eq. ([7]), with 8^ b replaced by A^ b , 
and the coefficients depend on the different form-factor normalizations Ri(l) (i = 0, 1, 2, 3) and 
the slope p 2 , see Appendix [D] for details. Here, again, we use inputs extracted from the decays 
involving light leptons where possible, while for the remaining form-factor normalization i?3(l) 
we adopt a value calculated in the framework of Heavy Quark Effective Theory (HQET) [33] . 

2.2 Overview of the A2HDM 

The 2HDMs extend the SM with a second scalar doublet of hypercharge Y — |. Thus, in 
addition to the three Goldstone bosons, they contain five physical scalars: two charged fields 
H ± and three neutral ones (p® = {h, H, A}. The most generic Yukawa Lagrangian with the SM 
fermionic content gives rise to tree-level FCNCs, because the Yukawa couplings of the two scalar 
doublets to fermions cannot be simultaneously diagonalized in flavour space. The non-diagonal 
neutral couplings can be eliminated by requiring the alignment in flavour space of the Yukawa 
matrices 21 ; i.e., the two Yukawa matrices coupling to a given type of right-handed fermions 
are assumed to be proportional to each other and can, therefore, be diagonalized simultaneously. 
The three proportionality parameters are arbitrary complex numbers and introduce new sources 
of CP violation. 

In terms of the fermion mass-eigenstate fields, the Yukawa interactions of the charged scalar 
in the A2HDM read |2l 



/2 

C% = ~ H + {u [q d VM d V R - ? u M U VV L ] d + Q vM{P R l} + h.c. , (9) 

where Vr : l = are the chirality projectors, and (/ = u, d, I) the proportionality param- 
eters in the so-called "Higgs basis" where only one scalar doublet acquires a non-zero vacuum 



expectation value. The CKM quark mixing matrix V 13 remains the only source of flavour- 
changing interactions. All possible freedom allowed by the alignment conditions is determined 
by the three family- universal complex parameters which provide new sources of CP violation 
without tree-level FCNCs |21|. Comparing Eqs. ^ and ([2]), one obtains the following relations 
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between the A2HDM and the general scalar NP parameters: 



9l 



™ M 2 H± ' 



9T dl 



-Wi 



* m g d m l 



M 2 H± 



(10) 



The usual models with natural flavour conservation (NFC), based on discrete Z2 symmetries, are 
recovered for particular (real) values of the couplings ^y; especially qi = q = — l/q u = — tan/3 
and <; u = qi = q = cot (3 correspond to the Type-II and Type-I models, respectively. 

Limits on the charged-scalar mass from flavour observables and direct searches depend 
strongly on the assumed Yukawa structure. The latest bound on the Type-II 2HDM charged 
Higgs from B X s7 gives M+ > 380 GeV at 95% confidence level (CL) pi. Within the 



A2HDM on the other hand it is still possible to have a light charged Higgs |24,25 . Assuming 



that the charged scalar H + only decays into fermions Uidj and l + vi, LEP established the limit 
M H ± > 78.6 GeV (95% CL) |35|, which is independent of the Yukawa structure. A charged 



Higgs produced via top-quark decays has also been searched for at the Tevatron 36 , 37 and 



the LHC 38,39 ; these searches are, however, not readily translatable into constraints for the 



model parameters considered here. It should be noted that the charged-scalar mass enters only 
in combination with the other couplings and, therefore, its size does not affect directly our 
results at this level. 



3 Results and discussions 

In Table [I] we summarize our predictions within the SM for the various semileptonic and leptonic 
decays considered in this work, using the hadronic inputs quoted in Table [3] (parameters that 
do not appear in this table are taken from [ll]). The rates for leptonic D, K and n decays 
are obtained from Eq. ^ with appropriate replacements, while the ratio of r — > K/nu T decay 
widths is given by [24] 

r(r -» Ku) f l-m 2 K /m 2 \ 2 
T(r — > nv) \ 1 — m 2 /m 2 J 

One can see that, apart from R(D) and R(D*), all the observables are in agreement with their 
SM predictions. While for the decays involving only Di a \ , K, and tt mesons no large effect could 
be expected because of the relatively small quark masses involved, this is equally true for the 
influence of gf? T . Contrary to that expectation, however, the data on R(D) and R(D*) indicate 
a large value for this coupling. This not only renders especially models with NFC incompatible 
with the data, but also poses a problem in more general models. 

We start by analyzing the constraints on the A2HDM parameters from the decays listed in 
Table [Tjpj In contrast to the models with NFC, the observables involving 5-meson decays (R(D), 

3 We do not take into account the experimental correlation between the measured values of R(D) and R(D*) 
given in [l] . It is reduced to — 19% when averaging with the Belle data and does not affect our results significantly. 
Moreover, the BaBar fit is sensitive to the assumed kinematical distribution, which is modified by the scalar 
contribution. While BaBar has already performed an explicit analysis within the Type-II 2HDM, it would be 
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Table 1: Predictions within the SM for the various semileptonic and leptonic decays discussed in 
this work, together with their corresponding experimental values. The first uncertainty given always 
corresponds to the statistical uncertainty, and the second, when present, to the theoretical one. 



Observable SM Prediction Exp. Value Comment 



R(D) 






0.296±^ ± 0.015 


0.438 ±0.056 
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0.354 ±0.026 




our averaj 
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(5.54 ±0.24) x 10" 
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R(D*) and Br (B — > tv)) can be consistently explained in the A2HDM. However, the resulting 
parameter region excludes the one selected by the leptonic -D( s )-meson decays. More generally 
speaking, models fulfilling the relations 

(a) gl^/gf^' = m qu m l /{m ql m v ) and (b) g^ 1 / ' gf^' = m qd rm/ \m q ,m v ) (12) 

are in conflict with the data. Removing R(D*) from the fit leads to a consistent picture in the 
A2HDM; in this case, however, the SM is also globally consistent with the data, as the tension 
in R(D) is "distributed" over the remaining observables. Models with NFC remain disfavoured 
compared to the SM. These observations lead us to consider the following scenarios: 

• Scenario 1 (Sc.l) is a model-independent approach where all couplings g^R 1 are assumed 
to be independent. One possible realization is the 2HDM of Type III. This implies that 
the effective couplings 5 l cb and A l cb in the two semileptonic processes can be regarded 
as independent. Therefore, predictions for the additional observables in B — > D(D*)tu 
follow in this case only from R(D) (R(D*)). 



• In scenario 2 we assume the relations in Eq. (12) to hold for qa = b, while processes 
involving only the first two generations are regarded as independent. When considering 
only the constraints from R(D), R(D*) and Br(B — > tv), the couplings in the A2HDM 
fulfill this condition; we will assume this form in the following for definiteness. 

• For scenario 3 we discard the measurement of R(D*) as being due to a statistical fluc- 
tuation and/or an underestimated systematic effect, leaving us with a viable A2HDM. 
From a global fit to all the other measurements we then obtain predictions for the new 
observables in B — > D^rv as well as R(D*). 

useful to analyze the experimental data in terms of the more general complex parameters A l cb and S l cb , to make 
the inclusion of this effect possible in the future. This modification is, however, only relevant for large values of 
the scalar couplings, which are excluded in the scenarios 2 and 3 discussed below. 
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The ratios in Eq. (12) might of course also be changed for I — r or I = /i. However, because 



of the smallness of m M , this would not be visible in any of the observables considered here. 
A way to test this option is to consider the ratio Br(B~ — > t~v t )/By(B~ — > fi^u^), which is 



independent of NP in the scalar sector if Eq. ( 12 ) is fulfilled 



In Fig. [T]we first show the allowed regions in the R(D)-R(D*) plane for the different scenar- 
ios. Scenario 1 is just reflecting the experimental information, while the additional constraint 
from B — > tv already excludes part of that area in scenario 2. For the third scenario, the 
tension with the measurement of R(D*) is clearly visible; the allowed range includes the SM 
range and values even further away from the measurement, thereby predicting this effect to 
vanish completely in the future if this scenario is realized. The value implied by the fit reads 
R(D*) Sc _ 3 = 0.241 ± 0.003 ± 0.007. 



of 




Sc.3 



0.4 0.5 

R(D) 




< o 




Figure 1: Allowed regions in the R(D)-R(D*) (left), complex S l cb (center) and A l cb (right) planes at 
95% CL, corresponding to the three different scenarios. See text for details. 

In the second and third plot in this figure we show the corresponding allowed parameter 
regions in the complex 5 l cb (left) and A l cb (right) planes at 95% CL. Here the strong influence of 
the leptonic decays becomes visible again, excluding most of the parameter space of R(D) in 
the 5 l cb plane, and driving the fit far away from the region indicated by the R(D*) measurement 
in the A l cb plane. To examine this effect further, we plot the individual constraints in the 
^i^u,d/ / M] I ± planes in Fig. [2j Here it is seen explicitly that the conflict lies mainly between 
the leptonic charm decays and R(D*). It is also noted that, in order to accommodate the 
data in scenarios 1 and 2, a large value for |w?f|/M|- ± ~ O(10 _1 ) GeV -2 is needed. This is 
not a direct problem in these scenarios, as the different couplings are not related. It would 
however point to a very strong hierarchy between the charm and the top couplings or a very 
large value of the leptonic coupling, as the constraints from observables involving loops like 
Br(Z — > bb) imply \q u \ < 1 for the top coupling 



24 



In the A2HDM, the combined constraints 
from leptonic r decays and these processes require a small value for \q u q*\/M^ I± : under the 
assumptions that |<^| < 50 and the charged-scalar effects dominate the NP contributions to 
Z — > bb, \<; u q*\/M^j± < 0.005 GeV~ 2 was obtained in [24]. Similar bounds also arise from 
considering the CP- violating parameter ex in K°-K° mixing and the mass difference Am^o 
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in B°-B° mixing 24 . This is however compatible with our results above, see Fig. [2J As the 
focus in this article lies on tree-level contributions, and loop induced quantities have a higher 
UV sensitivity, we refrain from including these constraints explicitly here. 
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Figure 2: Constraints in the complex <^£* /ML (left) and ? M ^*/M^ ± (right) planes, in units of 
GeV~ 2 , from the various semileptonic and leptonic decays. Allowed regions are shown at 95% CL for 
different combinations of the observables. 



Having in mind the above scenarios, our main concern is whether they can be differentiated 
with forthcoming data. In addition, the basic assumption of only scalar NP contributions in 
these decays can be questioned. We identify several combinations of observables which will 
signal the presence of additional contributions. In the following section we discuss how future 
measurements of additional observables in B — > D^rv decays, especially their differential 
distributions, will provide useful information to address these questions. 



4 Observables sensitive to scalar contributions 

We now proceed to analyze the additional observables in B — > D^*>tu t decays that provide an 
enhanced sensitivity to scalar NP contributions. Apart from the differential rates, these are 
the forward-backward and r-spin asymmetries in the considered decays, and for B — > D*tu 
additionally the longitudinal polarization fraction of the D*. These observables have been 
considered in the past by various authors [5^[6] [l^[M|[32^[42}]46] , addressing their sensitivity to 
NP contributions. 

Experimentally, while one expects more information on B — > D^*'t~u t decays in the near 
future from Belle, BaBar, and also LHCb, most of these observables will be accessible only at a 



Super-Flavour factory (SFF) (47 , 48 , as their study requires more statistics than the branching 
ratios and the inclusion of the information from the correlated B meson. The precise sensitivity 
of future experiments to the different observables has, however, not yet been determined. 
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Table 2: Predictions for the q -integrated observables both within the SM and in the different scenar- 
ios. The observables have been integrated from q 2 min = m 2 to q 2 max = {tub — m D (*)) 2 . The first error 
given corresponds to the statistical uncertainty, and the second, when given, to the theoretical one. 

Observable SM Prediction Scenario 1 Scenario 2 Scenario 3 

R L (D*) 0.115 ±0.001 ±0.003 0.217 ±0.026 0.223±g;g^| ± 0.006 0.1041°;°°^ ± 0.003 

A\(D) -0.304 ±0.001 ±0.035 — 0.55±g;J2 -°- 55 -ao4 -0.55+g;^ ± 0.01 

A\(D*) 0.502±g:gg| ± 0.017 0.06±g;Je - 04 iao3 ± °- 01 0.57±g$ ± 0.02 

Ag(D) 0.3602±g;ggg| ± 0.0022 0.03±g;gJ db 0.30 -0.21+°;^ ± 0.06 0.36±g;gJ + 

A e (D*) -0.066 ±0.006 ±0.009 — 0.136±g;ggj| ± 0.222 0.081±g;gg| ± 0.009 — 0.146±g;gfg =b 0.021 

1 Note that the lower tail is rather long, due to a suppressed local maximum. 



For simplicity, we discuss below these additional observables without considering the sub- 
sequent decays of the final r and D*. While studying the q 2 spectra of the observables has 
the advantage in identifying potential NP contributions and their Dirac structure, this gener- 
ally requires very high statistics, which might not be available even in the early stages of a 
SFF. Therefore, we shall analyze both the q 2 spectra and the g 2 -integrated observables. Our 
predictions for the latter are given in Table [2j both within the SM and in the three different 
scenarios defined before. Note that we do not consider isospin breaking; the observables shown 
are always isospin-averaged. Note furthermore that model-independent analyses similar to our 



scenario 1 have been performed in Refs. 5l 6 , 16 



4.1 The differential decay rates 



First, we obtain the singly differential rates by summing in Eqs. (36) and (43) over the r 
±1/2, and performing the integration over cos#: 
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(13) 



m T 
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(14) 



Again, normalizing to the decays with light leptons reduces the theoretical error: 

_ dB,{B -> D^r-v T )/dq 2 
H °^ q >- d Br(B^D(*)£-ve)/dq 2 ' U5j 

Note that in order to obtain the expression for R D (*) from this, numerator and denominator 
have to be integrated separately. This should be kept in mind as well for the other quantities. 

Our predictions for these observables, within the SM and in the three different scenarios, 
are shown in Fig. pi For Rf)(q 2 ), we show in addition the binned distribution in five equidistant 
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Figure 3: The q 2 dependence of the differential branching ratios (left) and R D (*)(q 2 ) (right), both 
within the SM (grey) as well as in scenario 1 (red), scenario 2 (orange), and scenario 3 (yellow). The 
binned distribution for Rn(q 2 ) is also shown. 



q 2 bins, as the ratio does diverge at the endpoint, where however both rates vanish. From these 
plots, we make the following observations: 

• As expected, the uncertainty due to the hadronic form factors (the grey shaded band) in 
R D (*)(q 2 ) is significantly reduced compared to that in the differential branching ratio. 

• In scenario 2, relatively large deviations from the SM are predicted for almost the full 
range in Ro*{q 2 ) and for high q 2 in Ro(q 2 ). The predicted q 2 spectra in scenario 3 are 
close to those of the SM, especially for the D* decay mode. However, there is still room 
for differences, especially in Ru(q 2 ), and furthermore in this scenario the distribution in 
RD*(q 2 ) lies preferably below the SM one, which should differentiate this scenario clearly 
from the first two. 

• Compared to the SM prediction, the peak of the differential branching ratio in the NP 
scenarios (especially in scenario 2) is shifted to higher and lower values in q 2 for the D and 
the D* decay mode, respectively. This is characteristic of scalar NP contributions and 
should allow for a separation from models with different Dirac structure. The reason for 



that is the following: while both of them are explicitly proportional to q (see Eqs. (32) 



and (39)), the charged-scalar contribution to Rd*(q ) is in addition proportional to the D* 
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momentum |p|, which vanishes at the endpoint g^ ax = (rag-mfl.) 2 , rendering its relative 
contribution maximal for intermediate values of q 2 , while the one to Rn{q 2 ) continuously 
increases with q 2 . The relative suppression of the terms proportional to the r mass by 
the D* momentum furthermore renders R,D*{q 2 ) finite everywhere, while Ro{q 2 ) diverges 
at the endpoint However, as both the rates for the r and the light lepton modes 
vanish there, this does not influence the experimental extraction: when calculating the 
contributions to different bins, all integrals remain finite. These characteristic features 
are illustrated by the right two plots. 

As the charged scalar only contributes to the helicity amplitude Hot in B — > D*tu t decays, 
an increased sensitivity is expected by studying the case with a longitudinally polarized D* 
meson in the final state, where the transverse helicity amplitudes are no longer relevant. For 
this purpose, we define a singly differential longitudinal decay rate |5]: 



dY L T G 2 F \V cb \\p\q 2 / m 



dq 2 967r 3 m 2 3 \ q 



,2\ 2 



,2 



Iff |2 / 1 , m A , 3m r \rr ,2 

\H 00 \ (1+2^) +2?\ H ot\ 



(16) 



In analogy to Rp* (g 2 ), it is again advantageous to consider the ratio with the r mode normalized 
to the light lepton mode: 

p* (i 2,_ dV L T /dq 2 

Rdq } - drfm ■ ( } 

It is important, however, to note that within our NP framework this is not an independent 
observable. As long as we consider only additional scalar operators, the difference 

X.iq 2 ) = R D ,(q 2 ) - Rl(q 2 ) (18) 

is independent of NP effects. A measurement of this observable serves, therefore, as a cross- 
check for the effect in Rp* and gives us information on whether scalar NP operators are sufficient 
to describe the data. This observation is reflected in Table [2] and in Fig. |4j where we show the 
predictions for Rl(D*) and R* L (q 2 ) both within the SM and in the three different scenarios; the 
results are analogous to the ones for R(D*) and Ro*{q 2 ) discussed above, but clearly exhibit 
an increased sensitivity to the scalar NP effect. 



4.2 The r spin asymmetry 

Information on the r spin in semileptonic i?-meson decays can be inferred from its distinctive 
decay patterns (5j[l6j|42j|46] . Therefore, we consider here the r spin asymmetry defined in the 
t-v t center-of-mass frame: 

DW( 2 = dT DM [X T = -1/2]/ 'dq 2 - dT DM [X T = +1/2]/ 'dq 2 

X {<1 ' dT DM [X T = -l/2]/dq 2 + dT DM [Xr = +l/2]/dq 2 ' 1 ' 



4 In fact, this behavior is an artifact of setting mi = 0. The actual value is ~ m\jm\ 
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Figure 4: Predictions for R* L (q 2 ) both within the SM and in the three different scenarios. The other 
captions are the same as in Fig. [3| 



where the polarized differential decay rates are obtained after integration over cos# of the 



doubly-differential ones given by Eqs. (36) and (43). Using the formulae presented in appen- 



dices [C] and [Dj we obtain explicitly 



2q*, 



3ml I tt 



\Ho\ 2 (l + 



00 1 



Ha. 



HP 



(1 



2?' 



3m| 
2g 2 



\H 



Of | 



00 I 



\H 



++\ 



+ \H_ 



'1 _l_ m ? ^ _1_ 3m ? I fT" 12 



(20) 



Again, these two observables have the same dependence on scalar NP contributions as the 
differential rates; this observation follows from the combinations 



X? (q 2 ) = R D {q 2 ) (A»(q 2 ) + 1) and (q 2 ) = R D *{q 2 ) (Af(q 2 ) + 1 



D* 



(21) 



being independent of 5^ b and A^ b , respectively. However, because of the different normalization 
and systematics in this case, a future measurement would give important information on the 
size and nature of NP in B — > D^*'tu decays: like for Xi(g 2 ), any deviation from the SM value 
of these combinations would indicate non-scalar NP. Our predictions for the two asymmetries 
when integrated over q 2 are given again in Table [2] the predicted ranges for the differential 
distributions are shown in Fig. M The correlation between A\(D*) and Rl(D*), following from 



Eqs. (18) and (21), is furthermore illustrated in Fig. 6 where we show the predicted values 



for the two observables for the SM and the three scenarios, yielding in every case a very small 
band, corresponding to the hadronic uncertainties. 

The integrated asymmetries span a rather large range. For A\(D), due to the common 
input, a differentiation between the different scenarios seems very difficult, while for A\(D*) 
at least the separation of the SM and scenario 3 on the one hand and scenarios 1 and 2 on 
the other hand is very clear. Importantly, the predictions for A\(D) in all scenarios are clearly 
negative, providing another option to potentially exclude the SM and only scalar NP at the 
same time. Also largely negative values for A\(D*) are excluded in all scenarios. 
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Figure 5: The q 2 dependence of the r spin asymmetries A®(q 2 ) (left) and (q 2 ) (right). The 
other captions are the same as in Fig. [3| 



0.2 



Figure 6: Predictions for A\{D*) vs. Rl(D*) both within the SM (grey) and in the three scenar- 
ios ((l)-red, (2)-orange, (3)-yellow), from a global fit including the appropriate observables. 



The separation of the different models improves, once differential distributions are consid- 
ered: from Fig. [5] we observe very distinct patterns for the different scenarios, especially for 
scenario 2. This scenario will be clearly distinguishable from the SM and scenario 3, once the 
necessary experimental precision is reached. Scenario 1, on the other hand, has again possi- 
bly very large effects, but might be close to any of the other scenarios, including the SM. A 
characteristic feature of this observable is the zero-crossing point: for A®(q 2 ), it is absent for 
the SM and scenario 3, while it appears likely for scenario 2 to have one, and also scenario 1 
has that option. Within the SM, the observable A®* (q 2 ) crosses the zero at q 2 = 3.66 ± 0.04; 
compared to the SM case, the zero-crossing point occurs at significantly higher values of q 2 for 
scenario 2, while most likely at lower values of q 2 for scenario 3. This indicates that measuring 
the zero-crossing point of the r spin asymmetries can be a useful probe of the flavour structure 
of the charged scalar interaction. 
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4.3 The forward-backward asymmetries 

Finally, we discuss the forward-backward asymmetries defined as the relative difference between 
the partial decay rates where the angle 6 between the and r three-momenta in the r-v T 
center-of-mass frame is greater or smaller than 7r/2: 

nM „ f ° d cos 9 (d 2 T? M ldq 2 d cos 6) - Q d cos (d 2 T? M ldq 2 d cos 6) 

A DM ,2) = J=l V JJo V r__^ )_ , , 

9 Kq ' dT?™ dq 2 1 ' 



Using Eqs. (36) and (43), we arrive at the following explicit expressions: 

3m 2 Re(H H*) 



2<7 2 \H \ 2 (l + §) + ^\H t \ 



A o* , 2) = 3 |i/ ++ | 2 -|i/-| 2 + 2^Re(ffoo^) 

4(|F ++ |2 + | J ff__|2+| J ff 00 |2)(i + g) + ^ W - 



In terms of a model- independent determination of NP parameters, i. e. scenario 1, this is the key 
observable to determine A T cb and 5 T cb . The reason, as mentioned before, is that the observables 
R* L (q 2 ), R D (,){q 2 ) and A^ * (q 2 ) do not give independent information, see Eqs. (18) and (21). 



The forward-backward asymmetry Af * (q 2 ) is therefore the only independent constraint in the 
complex Sl b (A^ b ) plane. Our predictions for this observable are given in Tableland shown in 
Fig. [7j Furthermore, its correlation with the r spin asymmetry for the two modes is shown in 
the first two panels in Fig. [8} It is clearly seen that the correlation is much weaker in this case, 
especially in scenario 1, where the only influence stems from the restriction on and |A£ 6 |. 
However, the pattern of a more SM-like A2HDM prediction and strongly shifted predictions 
from scenarios 1 and 2 is repeated. Regarding the differential distributions, within the SM, 
the observable A®(q 2 ) does not cross zero, while this becomes possible for scenarios 1 and 2. 
A®* (q 2 ) has a zero-crossing point at q 2 = 5.67 ± 0.02 in the SM, for which again large shifts 
are possible with NP, and it might even vanish in scenario 3. Large deviations from the SM 
expectations, especially in scenario 2, are therefore still possible for this observable. 

In order to illustrate the impact of a possible future measurement of this observable, we 
exemplarily show in the right panel in Fig. [8] the resulting constraint in the A^ fe plane, together 
with the one from R(D*) as measured at the moment. The Ag(D*) constraints drawn in lighter 
colours correspond to an uncertainty of 10%. For the darker constraints, an improvement 
by a factor of 2 has been assumed compared to the lighter ones. Furthermore, an index 
'SM' indicates the measurement chosen to be compatible with the SM, while the index 'NP' 
corresponds to measurements excluding the SM, but compatible with scenario 1. As can be 
seen, such a measurement would allow to exclude a large part of the parameter space in the 
model-independent scenario 1, as well as constrain the other scenarios further. Furthermore, as 
mentioned before, the two constraints could also miss each other in that plane, indicating NP 
with a different Dirac structure. This possibility exists of course also for the other observables 
discussed above. 
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Figure 7: Theq 2 dependence of the forward-backward asymmetries \q 2 ) (left) and A®* (q 2 ) (right). 
The other captions are the same as in Fig. [3| 

Additional information on A£ b could obviously be obtained from a measurement of the 
B~ — > t~v t rate. With the NP influence being determined by A^ b , this rate is clearly predicted 
to be different from the SM in scenarios 1 and 2, while close to the SM in scenario 3. However, 
this mode is extremely hard to be measured experimentally 




A e (D) A e (D ) Re(A' cb ) 



Figure 8: Prediction for Ag(D^) vs. A\(D(*>) for the SM (grey), and the three scenarios ((l)-red, 
(2)-orange, (S)-yellow), from a global fit including all appropriate observables. The right plot shows 
the possible impact of future measurements on the complex A£ fe plane (see text). 



5 Summary 

In this paper, motivated by the recent experimental evidence for an excess of r-lepton produc- 
tion in exclusive semileptonic 5-meson decays, we have performed a detailed phenomenological 
analysis of b — )■ qr~v T (q = u,c) transitions within a framework with additional scalar opera- 
tors, assumed to be generated by the exchange of a charged scalar in the context of 2HDMs. 

While the usual Type-II scenario cannot accommodate the recent BaBar data on B — > 
D^t~u t decays, this is possible within more general models, in which the charged-scalar cou- 
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plings to up-type quarks are not as suppressed. An explicit example is given by the A2HDM, in 
which the B — > D^*'r~v T as well as the B~ — > t~v t data can be fitted. However, the resulting 
parameter ranges are in conflict with the constraints from leptonic charm decays, which could 
indicate a departure from the family universality of the Yukawa couplings ^/ (/ = u, d, I). 

These observations led us to define three scenarios for scalar NP, in which we incorporated 
information from R(D^) (Sc.l), B decays (Sc. 2), and all available data from leptonic and 
semileptonic decays apart from R(D*) (Sc. 3). We showed that these scenarios can be differ- 
entiated by coming data, using information e.g. from differential decay rates and/or spin and 
angular asymmetries. These observables therefore allow to verify this hint for NP in semilep- 
tonic decays, and gather additional information on its precise nature. Furthermore we pointed 
out several combinations of observables independent of this kind of NP, as well as common 
characteristics, which will allow additionally to test for the presence of NP with other Dirac 
structures. The coming experimental analyses for these modes will therefore be an important 
step in our quest for NP. 
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Appendix 



A Input parameters and statistical treatment 

Bounds on the parameter space are obtained using the statistical treatment based on frequen- 



tist statistics and Rfit for the theoretical uncertainties 49 , which has been implemented in the 



CKMfitter package |l2). To fix the values of the relevant CKM entries, we only use determi- 
nations that are not sensitive to the scalar NP contributions 12,50 52 . Explicitly, we use the 



V u ,i value extracted from super-allowed nuclear /3 decays and the CKM unitarity to determine 
V us = A. The values of \V u b\ and \V c b\ are determined from exclusive and inclusive b —> utvi and 
b — > civ i transitions, respectively. Relevant hadronic input parameters are collected in Table [3j 
while quark and meson masses as well as any other relevant parameters that do not appear in 



this table are taken from 11 . 
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Table 3: Input values for the hadronic parameters, obtained as described in the text. The first error 
denotes the statistical uncertainty, and the second the systematic/theoretical. * This value includes the 



correction to the isospin limit usually assumed in lattice calculations 153, 661 



Parameter 



fB 3 /fB d 
fD s 

fD s /fD d 
Ik/ fir 

5, 



K 12/ ir 12 
cm 

tK2/ttt2 



A 

\V ub \ 

\Vcb\ 



Value 



(0.228 ±0.001 ±0.006) GeV 

1.198 ±0.009 ±0.025 
(0.249 ± 0.001 ± 0.004) GeV 
1.169 ±0.006 ±0.02 
1.1908 ± 0.0016 ±0.0104t 
-0.0069 ±0.0017 
0.0005 ± 0.0053 



0.2254 ±0.0010 
(3.51 ±0.11 ±0.02) x 10~ 3 
(40.9 ± 1.1) x 10~ 3 



Comment 



|54[j56| 



54 



55 



67 



57 



69 



54 


55 


57 


58 


54 


55 


57 


59 




59 


61 






62 


66 





70 



52 



52 



G{l)\V cb \ 

^\b^DIv 



^(1)1^*1 
fll(l) 

^2(1) 
#3(1) 
P 2 



1.186 ±0.036 ±0.041 
(42.64 ± 1.53) x 10~ 3 
0.46 ±0.02 



(35.90 ±0.45) x 10~ 3 
1.403 ±0.033 
0.854 ± 0.020 
0.97 ±0.10 
1.207 ±0.026 



52 



52 
|7||l7|p8' 



52 



52 



52 



33 



52 



The plots for the differential observables are obtained using the allowed NP parameter 
ranges from a different fit. As the latter already include uncertainties from the hadronic input 
parameters, we do not vary them again additionally. 

B Kinematics for semileptonic decays 

Within the SM, the squared matrix element for the decay B{pb) — > D^(p D (*), X D (*))l(ki, Xi)u(kp) 



can be written as 31 32] 

\M{B D^W)\ 2 = \(D^W\£ eS \B)\ 2 = L, U H^ , (24) 

where the leptonic (L^) and hadronic (H^) tensors are built from the respective tensor prod- 
ucts of the lepton and hadron currents. Using the completeness relation for the virtual W* 
polarization vectors e M (±,0,t), one can further express Eq. (24) as 



\M{B -»■ D w /P)| 2 = L{™,n)H{m',n')g mm ,g nn ,, (25) 

m,m' ,n,n' 

where g mm > = diag(±l, —1, —1, —1), L(m, n) = L^ v e^(m)e* u (n) and H{m, n) = H^ v e* (m)e„(n). 
The two quantities L(m, n) and H(m, n) are Lorentz invariant and can, therefore, be evaluated 
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in different reference frames. For convenience, the hadronic part H(m,n) is usually evaluated 
in the 5-meson rest frame with the z axis along the D^*' trajectory, and L(m,n) in the l-v 



center-of-mass frame (i.e. in the virtual W* rest frame) [3T 32 . 

In the .B-meson rest frame with the z axis along the trajectory, a suitable basis for the 
virtual W* polarization vectors e M (±,0, t) can be chosen as [31] 

e„(±) = i(0,±l,-i,0), e M (0) = -^=(|p|,0,0,-g ) 



, q 2 

^(t) = ^=(g ,0,0,-|p|), (26) 



where go = ( m B — m2 D {*) + Q 2 )/^ m B and |p| = A 1/l2 (m|,m^ w ,g 2 )/2mB are the energy and 
momentum of the virtual W*, with q 2 = (jps — Pdm) 2 being the momentum transfer squared, 
bounded at mf < q 2 < (tub — m D (,)) 2 , and A(a, b, c) = a 2 + b 2 + c 2 — 2(ab + be + ca). Similarly, 
a convenient basis for the D* polarization vectors is 

e a (±) = T -^(0,l,±i,0), = — (|p|,0,0,^), (27) 

V 2 mo* 

where E D (*) = (m B + rn 2 D{t) — g 2 )/2m^ is the energy in the S-meson rest frame. 

In the l-v center-of-mass frame, which can be obtained by a simple boost from the 5-meson 
rest frame, the lepton and antineutrino four-momenta are given, respectively, as 

ki = (Ei,pi sin 0,0, p/ cos 9) , k p = (p h -p t sm9, 0, -p t cos 9) , (28) 

where E\ = (q 2 + m 2 )/2^y~q 2 , p\ = (q 2 — m 2 )/2y(f, and 9 is the angle between the and / 



three-momenta in this frame. The virtual W* polarization vectors e M (±,0,t) reduce to 31 32 



e„(±) = -J= (0, ±1, -i, 0) , e M (0) = (0, 0, 0, -1) , 
e,(t) = -^=<z„ = (1,0,0,0). (29) 



With the above specified kinematics, the explicit expression for L^ V H^ V in terms of the q 2 
dependent helicity amplitudes can be found in Refs. (5,31 . Using the equations of motion, 



the hadronic and leptonic amplitudes of the scalar and the pseudoscalar current can be related 
to those of the vector and the axial- vector current, respectively. Therefore, the scalar NP 
contributions can be considered together with the spin-zero component (Xw* = t) of the virtual 
W* exchange. 
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C Formulae for B —> Dlv 

In the presence of NP of the form g, the non-zero hadronic matrix elements of the B — > D 
transition can be parametrized as 



(D(p D )\cYb\B(p B )) 



(D(p D )\cb\B(p B )) 



U(q 2 ) 



{Pb+PdT- 



m B — m 



v 



h{q 2 



m 



B 



{D(p D )\SYb\B(p B )) 



m 



B 



m 



D 



(30) 
(31) 

nib — m c mb — m c 

where fn q are the running quark masses and the two QCD form factors f+(q 2 ) and fo{q 2 ) encode 
the strong-interaction dynamics. Contracting the above matrix elements with the virtual W* 
polarization vectors (26) in the B- meson rest frame, we obtain the two non- vanishing helicity 
amplitudes [5 31 : 

2m B \\)\ 



H (q 2 ) 
Ht(q 2 ) 



U(q 2 



m, 



m 



D 



fo(q 2 



l + 6i 



cb 



(m B - m D ) 



(32) 



where S l cb , defined by Eq. (J6|, accounts for the contribution from the charged scalar. 

It is customary to relate the QCD form factors f+(q 2 ) and fo(q 2 ) to the quantities G\(w 
and A(w) in the HQET [33] 

G\(w) 



U(q 2 



R 



fo(q 2 



D 



(33) 



where R D {*) - 
is defined as 



2 y /m B m DM /(m B + m DM ), r 



w 



(1 + w) 1 + r 

2 1 — r 

m D (*) /m B , and the new kinematical variable w 



m\ + m 2 D(t) - q 2 



v B -v DM = — — , (34) 

2m B m D M 

with v B and v D (*) being the four-velocities of the B and mesons, respectively. We approx- 
imate the scalar density A(w) by a constant value A(w) = 0.46 ± 0.02 (7, 17,28 and G(w) is 
parametrized in terms of the normalization Gi(l) and the slope p\ as (72j 

G x (w) = Gi(l) [1 - 8pl z(w) + (51 pi - 10) z[w) 2 - (2h2p\ - 84) z(w) 3 ] , (35) 

with z(w) 

Equipped with the above information, the double differential decay rates for B — > Dlu, with 
/ in a given helicity state (A; = ±1/2), can be written as 



d 2 r g [Az = -1/2] 
dq 2 d cos 6 

^r g [A, = +i/2] 

dq 2 d cos 9 



G% 


V cb \ 2 q 2 






m 


l7i 3 m 2 B 


q 2 




V cb \ 2 q 2 




nr[ 


m 


l7r 3 m 2 B 


q 2 



\p\\H (q 2 )\ 2 sm 2 6 



iPl 



ml 



q 



\H (q 2 ) cos 6 -H t (q 



2\|2 



(36) 



from which the total decay rate and the various g 2 -dependent observables can be obtained via 
summation over A; and/or integration over cos#. Owing to its lepton-mass suppression, the 
A; = +1/2 helicity amplitude is only relevant for the r decay mode. 
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D Formulae for B — > D*lv 

For the B — > D* transition, the hadronic matrix elements of the vector and axial-vector currents 
are described by four QCD form factors V(q 2 ) and A i^(q 2 ) via 

(D*(p D *,e*)\c^b\B(p B )) = 2lV } q2) e^iTjfoh, 

m B + rnc* 

(D*(p D *,e*)\c W b\B(p B )) = 2m* D A (q 2 ) e —^-q^ + (m B + m D *)A 1 (q 2 ) (^-^q^ 



Q- 
e -q 



m B + m D * 



(Pb+Pd*] 



m 2 B — m 2 D * 



(37) 



from which one can show that, using the equations of motion, the B — > D* matrix element for 
the scalar current vanishes while the pseudoscalar one reduces to 



{D*(p D *,e*)\c l5 b\B(p B )) 



nib + m c 



2m 



D 



nib + m c 



(D*(p D *,e*)\c^ l5 b\B(p B )) 
A (q 2 )e*-q. 



(38) 



Contracting the above matrix elements with the W* and D* polarization vectors in Eqs. (26) 
and (27), we obtain the four non-vanishing helicity amplitudes [5,31 



H ±± (q 2 ) 
Hoo(q 2 ) 
H 0t (q 2 ) 



{m B + m D *) A^q 2 ) =F — |p| V{q 2 



2niD* v q 2 
2m B \p 



m B + m D * 



(m B - m D , - q ) (m B + m D ») A t (q 



m B + m D . 



A 2 (q 2 ) 



A (q 2 ) 1 - A 



cb 2 ( ' 



m 



B 



(39) 



where A^, defined by Eq. Q, accounts for the contribution from the charged scalar. 

In the heavy-quark limit for the b and c quarks, the four QCD form factors V(q 2 ) and 



Ao,i ; 2{q 2 ) are related to the universal HQET form factor Ha^w) via 72 



V(q 2 



A n (q 2 ) 



D* 



h Al (w) 
h Al {w) 



A 1 ( 



A, 



Rd* 
Rq(w] 



R„ 

Rd* h Al (w) 



R2H 
Rd* 



h Al {w) 



(40) 



21 



where the w dependence of h Al (w) and the three ratios Rq^^w) reads [72] 

h Al {w) = h Al (l) [1 - 8p 2 z(w) + (53p 2 - 15) z(wf - (231p 2 - 91) z(wf] , 

Ro(w) = Ro{l) -0.11(w- 1) + 0.01(w- l) 2 , 

i?i(w) = i?i(l) -0.12(iu- 1) + 0.05(w- l) 2 , 

J2a(i£7) = Ra(l) -0.11(w- 1) - 0.06(w - l) 2 . (41) 

The free parameters p 2 , -Ri(l) and i?2(l) are determined from the well-measured B — > D*tv 
decay distributions [52] (£ = e,p), whereas for the parameter Ro(l), that appears only in 
the helicity-suppressed amplitude H ot , we have to rely on the HQET prediction for the linear 



combination 33 



R 2 (i)(l-r)+r\Ro(l)(l + r)-2] . . 

R 3 (l) = -^-^ (1 -r)2 = °- 97±0 - 10 ' ( 42 ) 

which includes the leading-order perturbative (in a s ) and power (l/m& )C ) corrections to the 
heavy-quark limit, plus a conservative 10% uncertainty to account for higher-order contribu- 
tions [5]. 

Finally, the double differential decay rates for B — > D*W, with / in a given helicity state (A/ = 
±1/2), can be written as 



rf 2 r^[A, = -l/2] = G 2 F \V cb \ 2 \p\q 2 / 

dq 2 d cos 9 25Qir 3 m 2 B \ q 2 



B 

x [(1 - cos6») 2 |i^ ++ | 2 + (1 + cos^) 2 |i^r I 2 + 2sin 2 6' |i^oo| 2 ] , 



d 2 r^[A, = +l/2] G 2 F \V cb \ 2 \p\q 2 ' 



1 



,2 



dq 2 d cos 9 256n 3 m 2 B \ q 2 J q 

x [sin 2 9 (\H ++ \ 2 + \H__\ 2 ) + 2 \H ot -H 00 cos 9\ 2 } , (43) 

which are the starting point for the total decay rate, as well as the additional observables 
considered. 
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